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O ' Abstract 

o : 

\ The a resonance was observed as a conspicuous tt+tt^ peak in hadronic decays 

^ ■ like J/if) ^ TX^TX^ uj or 7r+7r^7r+. The phase of the a — * tt+tt" ampli- 

■ tude, extracted from production data within the conventional isobar model, is 

^ . assumed to coincide with that in tttt elastic scattering. We check the validity 

! of this assumption by using Lehmann-Symanzik-Zimmermann (LSZ) reduc- 

CN \ tion and unitarity. The rescattering effects in the final three-particle states 

^ ■ are shown to generate a correction to the phase given by a naive application 

lO ■ of Watson theorem. We briefly discuss the implications of this result for the 

^ . pole determination from production data. 

§ : PACS: 13.25. Gv; 13.20.Fc; 13.75.Lb; ll.80.Et 

O ; 

Ph! 1 Introduction 

^ ■ The lowest scalar resonance a (or /o(600)) appears as a pole on the second 

Riemann sheet of the / = / = partial wave amplitude of vrvr elastic scattering 

! (we denote this wave as to(s)). Although a typical resonant behaviour is not 

^ \ seen, because the pole is far from the real axis and is compensated in the 

' physical region by the Adler zero, many determinations of the sigma pole 

are based on vrvr scattering [T]. However, the pole was usually extracted 
from parametrizations valid along the physical region. The predictions are 
therefore affected by the large uncertainties of the analytic extrapolation to a 
distant point. Recently, a model-independent extrapolation into the complex 
plane, based on the Roy equation for to(s), led to a precise prediction of the 
pole position [2]. 

The a resonance was also seen as a peak in BES II data on J/ip vr+vr^u; 
[S] and in the data on —>■ vr"''vr~vr"'" reported by E791 Collaboration 
The conspicuous sign in production processes is explained by the absence of 
the Adler zero El- It is of interest to compare the pole determinations 
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from production processes and tttt elastic scattering. In the present paper we 
consider some issues related to this problem. 

To illustrate the discussion we consider the strong decay 

J/^(p) ^ 7r+(pi) + n'{p2) + uj{p^) , (1) 

but our arguments apply also to the decay 7r+7r^7r+, and more generally 

to h ^ 71^71^ hi, where h and hi are hadrons. We define the Mandelstam 
variables 

S=iPl+ P2)\ t={pi+ p^)\ U={p2+ Pzf. (2) 

which satisfy s + t + u = mjy^ + 2m^ + m5. The physical region of the process 
(P) corresponds to s > 4m^, t > (m^ + m^)'^ and u > (m^ + m^^)^. Since 
some particles have nonzero spins, a decomposition in Lorentz covariants 
is required. For simplicity, we neglect this complication and consider an 
invariant amplitude A{s,t) as a function of the Mandelstam variables (0). 

In the conventional isobar model, the amplitude of the decay (P) is ex- 
pressed as a sum of isobaric resonances in various channels^. In a diagram- 
matic language, the three-body decay is assumed to be described by tree 
diagrams where the production of two-body final states proceeds via inter- 
mediate resonances. More exactly, the amplitude is written as jSl E] 

A{s, t) = Asis, t) + [Atit, s) + it^ u)] , (3) 

where the s-channel and t (M)-channel amplitudes are expanded as 

A,{s,t) = J2ai{s)Pi{cose,) , A{t,s) = Y.bi{t)Pi{coset) . (4) 

In these relations 6s{= 613) is the angle between the three-momenta of vr(pi) 
and uj^ps) in the rest system of the two pions, and Ot{= 612) the angle between 
the three-momenta of vr(pi) and vr(p2) in the rest system of T{{pi)uj. The 
lowest partial waves in ^ are assumed to be dominated by resonances. For 
the process (H)) the s-channel resonances a, /o(980) and /2(1270) contribute to 
the partial waves and 02, respectively, and 65'"(1235) appears in both S and 
D-waves of the t (u) channels. Keeping for simplicity only the contributions 
of a and bi and assuming Breit-Wigner parametrizations, one writes [3] 



ao{s) = — — pY^, bo{t) = — -i — : — . (5) 



complex constant accounting for the direct nonresonant interaction is sometimes 
added to the resonances (see 01 E|)- 
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In Refs. OEl, ao(s) is denoted as the a ^ nn amplitude. 

The phases A^- and A^^ appearing in (0) account for the interactions of 
au and irbi, respectively. In the conventional isobar model these phases 
are assumed to be independent on the Mandelstam variables. Moreover, by 
invoking Watson theorem jTj, the phase of the a rni amplitude ao(s) was 
assumed El to coincide, up to a constant, with the pion-pion phase shift 
6q appearing in the expression of the / = / = partial wave: 

t'ois) = ^Ms)e'<'^''> - 1} , (6) 



where p{s) = yl — Am^/ s. An equivalent formulation is to assume jE] that 
the denominator of the function ao(s) given in (jSJ coincides with the function 
D[s), appearing in the N/D method [H] for calculating the amplitude t^{s). 

The purpose of this letter is to check the validity of Watson theorem in the 
isobar model for decay processes. We recall that the theoretical difficulties of 
the three-particle decays are known since a long time. Anomalous singularities 
generated by rescattering effects and three-body dispersion relations were 
considered by several authors (see EDI and older references quoted therein) . 
In the present paper we investigate the phases of the amplitudes defined in 
the isobar model, using an approach based on LSZ reduction and unitarity. 



2 LSZ reduction and unitarity 

We start from the S'-matrix element of the process (0) 

Sft = (7r(pi) 7r(p2)^(p3); out \J/ip{p); in) . (7) 
After the LSZ reduction fTH, IT^ of the u meson we obtain 

Sfi = 6f, + — ^ / dxe'P'-^{7r{p,) Tiip^); out IvMlJ/^ip); in) , (8) 

where ps^ is the time component of ps and rj^^ = ICx4>uj{x) denotes the source 
operator {JCx is the Klein-Gordon operator and (p^j the interpolating field of 
the omega meson). In what follows we do not need the explicit expressions of 
the sources, but only the significance of the matrix elements involving them. 
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Using translational invariance rj^^ix) = e^^'^'r]i_j{0)e ^^'^ where P denote 
the momentum operator, we write (jHI) as 



Sfi = 6f,+ (27r)^(5 (pi +P2+P3-P) (7r(pi) n{p2) ; out |?7^(0) \J/ij{p);m) . 

(9) 

From the general expression of the S'-matrix in terms of the invariant amph- 
tude A{s,t), it follows that 

A{s,t) = -^={n{p^)n{p2);ont\T]u.mJ/ij{p)m) . (10) 

where JV = 4^j9o]5i,oP2,oP3,o is a normalization factor. In the same way we ex- 
press the invariant amplitude T(s, t) of the elastic scattering n^ki) + Tr{k2) — >■ 
7r(pi) + n{p2) as 

T{s,t') = ^=(vr(pi) |r/^,(0)|7r(fci)7r(fc2);in) , (11) 



where J\f' = 4ypi 0^2,0^1,0^2,0, and the physical domain in the s channel is 
defined by s = (pi + ^2)^ = (^i + h)'^ > 4m^, t' = {pi - ki^ < 0. 

By applying once more the LSZ reduction to the matrix element (jlUj) . we 
obtain: 

A{s,t) = , / dxe*^--^(xo)(7r(pOlk.(x),r/^(0)]|J/V^(p);in) , (12) 



^2p3,o2p2,l 

where riT,^{x) is the source of the final pion Tc{p2). 

As it is known, the LSZ formalism allows the analytic continuation of the 
amplitude A{s,t) in the complex planes of the Mandelstam variables, where 
the expression (jl2p defines a holomorphic function (an important ingredient in 
the proof is causality, i. e. the fact that the retarded commutator vanishes for 
spacelike values x"^ < 0). In what follows we only use the LSZ representation 
to derive the unitarity relation and explore its consequences. 

By inserting a complete set of states \n) in the two terms of the retarded 
commutator, the matrix element appearing in (|12j) writes as 

E [{7vipi)\Vn,ix)\n){n\U0)\J/i^) - {7rip,)\U0)\n){n\v.,{x)\J/i^)] ■ (13) 

n 

In the two particle approximation, the lowest states which contribute in the 
first sum are \n) = 1 71(^1)71(^2)), where the two pions have / = 0, while 
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in the second sum \n) = \7T{ki)uj{k2)) ■ The sum over intermediate states 
involves integrations upon the momenta of the on-shell particles and sums 
over polarizations. After imposing the translation invariance, we obtain from 



A{s,t) =iAf dx6{xo) 



^•^1 dk2 ip2-x+ipi-x—iki-x—ik2-x 

(27r)3 (27r)3 

K 

J J (27r)'^ [2'nY 

(7r(pi)|r7^(0)|7r(A:i)cu(fc2))(7r(A:i)cu(A:2)|r/.,(0)|J/V;(p)) 



(14) 



The presence of 9{xo) in the integral of Eq. ()14|1 leads to discontinuities of the 
amplitude across the real axis of the Mandelstam variables fUE]. According 
to the general prescription, the discontinuity is obtained formally from the 
expression ()14p through the replacement of i9{xQ) by 1/2 T2]. Then the 
integral upon x gives (27r)'^5(pi+p2 — ^i — ^2) in the first term, and (27r)'^(5(p — 
P2 — ki — /C2) in the second. Moreover, in the first term of (jl4p we recognize 
the product of the amplitudes of the processes n^ki) + vr(A;2) t^{Pi) + t^{P2) 
and J/i'ip) — > 7r(A;i) + 7r(A;2) + 07(^3), while in the second term appears the 
product of the amplitudes of the processes 7r(fci) +uj{k2) +^^^(^3) and 

J /Hp) 7r(p2) + vr(fci) + uj{k2). 

From this discussion, it follows that the amplitude A{s, t) has branch cuts 
for s > Anil and t > {171^, + m^Y . By reducing the pion 7r(pi) instead of n{p2) 
one obtains the branch cut for u > (m^r + mt^)^. Hence, in the physical region 
of the decay all the variables s, t and u are above the unitarity thresholds. 
Singularities in both the s and t variables are expected to occur in each of 
the two terms in (jl4j) . We recall that in the isobar model defined in ©-(01), 
the term As{s,t) has singularities only in the s variable, being holomorphic 
with respect to t, while At{t, s) has singularities only in t, being regular with 
respect to s. This shows the limitation of the isobar model. 

For definiteness we consider that the complete set of intermediate states 
are "out" states and recall that for one particle states the "in" and "out" sets 
are equivalent. Recalling the definitions (fTUI) and (fTTj) of the invariant ampli- 
tudes and focusing on the first term in (|T^. responsible for the discontinuity 
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with respect to the variable s, we obtain 

+ ,) _ A(s - >M)} = gij / T'(s, W) . (15) 

where P = Pi+P2 — ki — k2. The amphtudes are evaluated for s = (^1+^2)^ = 
{ki + /c2)^ and the momentum transfers t' = [pi — fci)^ and t" = {ki +^3)^, 
respectively. The integral (jl5j) is easily evaluated in the c.m.s. of the two 
pions pi + p2 = ki + k2 = 0. After the trivial integrations due to the delta 
functions, (fTSj) reduces to an integral upon the angular variables: 

L{A{s + te, t) - A{s -te,t)} = -^ J dn p{s)T* {s, t')A{s, t") , (16) 

where dfi = d0dcos^", 9" being the angle between the three momenta of 
T^{ki) and uj{p^) in the pion rest system. For the vrvr isoscalar amplitude 
T{s,t') we use the Legendre expansion 

00 

T{s,t') = 167r^(2/ + l)P,(cose')^|'(s), (17) 
1=0 

where 6' is the angle between the three momenta pi and ki. We stress that 
()16|) is a general unitarity relation, independent of the isobar model. Let 
us restrict now to this model, taking for the amplitude A{s,t") in (|16p the 
expression given in Using the well known relation jS] 

cos 9" = cos 9 cos 9' + cos sin 9 sin 9' , (18) 

the integral upon the angle is trivial. Recalling that only the first term in 
Q has a discontinuity for s > 4:ml, projecting onto the S'-wave and using 
the orthogonality of Legendre polynomials we obtain 

^{ao(s + ie) - ao(s - ie)} = p{s) (t[J(s))*ao(s) 

+^it'ois)r I dcos9"[Ait", s) + Au{u",s)]. (19) 
-1 

We recall that time reversal invariance implies the reality relation ao(s — ze) = 
a^ls + ie), from which it follows that the l.h.s. of (fT^ is real and equal to 
Imao(s) (if not otherwise specified, s is taken on the upper edge of the cut). 
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3 Watson theorem 

Neglecting the four-pion channel which opens very slowly, the elastic region 
extends up to the threshold for KK creation. Below this threshold viq{s) = 1, 
therefore the amplitude (jHl) becomes tg = e^^° sin 5q / p{s) . If we neglect the 
second term in the r.h.s. of (fTTHl we obtain 

^{ao(s + le) - ao{s - le)} = e'^^''^ sin5°(s) ao{s) . (20) 

This relation implies ao(s + ie) = ao{s — ie)e'^^^o(^) ^ which is equivalent to 
Watson theorem: the phase of aQ{s) is equal (modulo ivr) to the phase shift 
6q. Alternatively, writing [S] 

where N{s) has only a left hand cut for s < and D{s) a right hand cut for 
s > 4:m'^, a solution of ()20|) has the form: 

.„W = §|i|. (22) 

where the function C{s), real for s > Am^, is arbitrary. By the uniqueness of 
analytic continuation, this implies that t^{s) and the function a^i^s) have the 
same poles on the second sheet. We will come back on this point in the next 
section. 

If the second term in the r.h.s. of (fT^ is not neglected we obtain, instead 
of (j2ni), the more general relation 

Imao(s) = e-<^'^ sin5°(s)[ao(s) + h{s)] , (23) 

where 

h{s) = -j dcose"[At(t'\ s) + [Au{u'\ s)]. (24) 

-1 

We recall that the angles 6s and 6t in the expansions (jlj are expressed in 
terms of the Mandelstam variables, for instance: 



ml + ml + + ml^s-t 

cos Os = = , (25) 

2\p\Js/4-ml 



where \p\ = A"^/^(s, rrij^^, m^)/(2A/i) is the three momentum of J/ip (i^) in 
the rest system of the pions (here A (a, b, c) = a"^ + b"^ + — 2ab — 2ac — 26c). 
Using these relations and retaining in the expansion (0]) of + Au only the 
S-wave bo, parametrized as in (0), we have: 



iA 



^bi - ml - ml- iTb.rrib^ - ps^o^/s - 2\p\\k2\ 



(26) 

where 1^21 = — and p^^o = ^\p\'^ + ml, with \p\ defined below (j^ . 

If the isobar model contains also a nonresonant term jHIH], the function h{s) 
will include its contribution. 

From Eq. (j23|) we can calculate a correction to Watson theorem. To this 
end we impose time reversal invariance, which means that the r.h.s. of (j23j) 
must be real. By requiring that the imaginary part vanishes, we obtain 



Im 



smMs) - 5°(.)] = ' , , (27) 



|ao(s)| 

where $(s) is the phase of the production amplitude: 

ao{s) = \ao{s)\e''"~'\ (28) 

The relation ^I7\i gives a calculable correction to the phase predicted by the 
naive application of Watson theorem in the elastic region. Above the inelastic 
threshold s = Am\- the elasticity iil{s) in (jU)) drops very quickly below unity, 
and additional terms due to the KK intermediate states appear in the r.h.s. 
of (Uni). Since the unitarity sums for the scattering and the decay processes 
contain different contributions, the phase of ao(s) in the inelastic region may 
be quite different from 5q(s). 



4 Comments 

The above analysis shows that the phase of the a tttt amplitude cio('S), 
defined in the conventional isobar model for hadronic decays, is not exactly 
equal to the tttt phase shift, as one would think by a naive application of 
Watson theorem. In the isobar model, the complex constants multiplying the 
Breit-Wigner resonances describe the interaction of a resonance (a) with the 
third hadron {oj). The s-dependent correction $(s) — 5q(s) calculated above 
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is generated by the individual interactions with uj of each of the outgoing 
pions. Actually, the rescattering effect discussed above can be visualized by a 
triangular diagram, given in Fig. 1 of Ref . . As shown in jU] , this diagram 
is responsible for the appearance of anomalous singularities. In the present 
work we emphasized the influence of the rescattering effects on the phase of 
the 0" ^ TTTT amplitude defined within the isobar model. 

We notice that, if the total amplitude t") appearing in the r.h.s. of the 
unitarity relation ()16|) could be expanded in a series of Legendre polynomials 
Pi(cos^s), the standard evaluation of the integral |H] would lead to Watson 
theorem for each partial wave a;(s). In the case of the elastic vrvr scattering 
(or in decays like K n-nlv^ such an expansion is legitimate, since in the 
physical region of the s-channel the amplitude is a holomorphic function of 
t. On the other hand, for three-body decays like (0) a similar expansion is 
not possible, since Pz(cos6's), which are polynomials of t, fail to reproduce 
the branch cut along t > {m-,, + m^)^. The isobar model attempts to take 
into account the singularities in all channels, but, as we discussed above, it is 
too simplistic. As a consequence, the phase of the a ^ tttt amplitude ao(s) 
defined in this model is not exactly equal to the phase-shift 5q{s). 

From it follows that the magnitude of the phase difference $(5) — 5o(s) 
depends on the values of the parameters of the isobar model (the ratio Cb^ / 
and the difference A^,^ — Aq.). Since an overall constant phase is irrelevant, 
what really matters is the variation with s. The difference $(s) — 5g(s) 
might be smaller than the experimental errors^. However, it is important to 
emphasize that even a small phase difference may have an important influence 
on the pole determination. Indeed, an immediate consequence of our result 
is that the denominator of the function ao(s) in (0) should not be identical to 
the denominator D{s) appearing in the expression (j^lj) of t^{s). In Refs. jSlE] 
the a pole is extracted from a parametrization of the denominator of ao(s) 
in the physical region, supposed to be valid also in the complex plane. As 
mentioned in the Introduction, such a method is affected by the uncertainties 
of analytic continuation, which are large for a distant pole. If, in addition, the 
denominator of ao(s) and the denominator of the vrvr amplitude t^{s) differ 
in principle (even slightly) along the physical region, the position of the pole 
determined by the analytic continuation of the production data may be even 

^For the decay — > 7r+7r~7r+, where the statistics is rather low, fits of equal quality 
were obtained both with a phase of ao(s) very different from the tttt phase-shift j4j, and 
with a phase close to 5q 0. 
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more distorted. The effect discussed above migfit play a role in understanding 
the difference between the mass and width of the lowest scalar resonance a 
extracted from the BES data for J /ip ^ n^n'uj decay 13] , and the values 
derived from Roy equation for the / = Z = elastic vrvr amplitude [2]. 
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